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I.  INTRODUCTION 


RF  measurement  accuracy  is  a  fundamental  assessment  in  experiment  design, 
measurement  programs,  and  system  evaluations.  Projections  of  overall 
measurement  accuracy  are  typically  made  through  error  budget  assessments 
comprised  of  an  rss  sum  of  the  random  error  components  added  to  the  sum  of 
bias  error  components.  The  error  budget,  therefore,  requires  knowledge  of  the 
individual  error  mechanisms  and  quantification  of  their  first  and  second  order 
statistics.  Over  the  years,  significant  use  has  been  made  of  Gaussian 
statistics  for  error  projections  dating  to  the  original  application  by 
Woodward  (Ref.  1).  A  tutorial  presentation  of  this  analysis,  its  underlying 
validity,  and  applications  were  presented  by  Swerling  (Ref.  2),  and  the  extent 
of  such  applications  is  described  in  Ref.  3.  Further  extensions  of  this 
analysis  continue  into  the  present  time;  the  statistics  for  the  phase 
difference  between  two  vectors  have  recently  been  published  (Ref.  4).  These 
analyses,  however,  are  tied  to  the  fundamental  assumption  of  Gaussian 
statistics. 

When  errors  are  coherently  related  to  the  desired  signal,  the  assumption 
of  additive  Gaussian  error  statistics  conflicts  with  physical  interpretation . 
The  error  statistics  for  the  coherent  case  are  derived  based  on  the  familiar 
phasor  diagram  shown  in  Fig.  1,  where  the  true  value  is  set  to  unity  with  a  0° 
phase  angle  without  loss  of  generality,  and  the  error  has  a  constant  relative 
amplitude  "a"  and  phase  angle,  a,  which  is  assumed  to  be  equally  likely  and 
uniformly  distributed  from  0  to  360°.  The  physical  model  of  this  situation  is 
particularly  appropriate  with  today's  popularity  for  swept  frequency 
measurements.  Generally  the  magnitude  of  the  error  component  varies  slowly 
with  frequency  changes  while  the  phase  between  the  true  and  error  components 
rotates  with  frequency  changes.  This  phasor  picture  of  the  errors  has  had 
most  application  and  familiarity  in  terras  of  the  peak-to-peak  values  of  the 
error.  However,  the  first  and  second  order  statistics  for  this  distribution 
have  had  relatively  little  investigation  and  are  believed  to  be  original  in 
this  work.  The  simple  phasor  diagram  is  adequate  for  narrow  bandwidth 


signals;  for  broader  bandwidth  waveforms,  the  error  component  needs  to  be 
weighted  by  the  autocorrelation  of  the  waveform  with  the  appropriate  time 
delay. 


Fig.  1.  Phasor  Diagram  of  Coherent  Error  Analysis 

If  both  the  amplitude  and  phase  of  the  error  component  can  be  determined, 
then  the  error  can,  of  course,  be  quantified  deterministically  by 
calibration.  Indeed,  this  technique  is  well  developed  in  modern  network 
analyzer  practice,  where  vector  measurement  techniques  used  in  conjunction 
with  a  set  of  reference  standards  obtain  the  full  scattering  matrix 
characterization  of  the  measurement  system  (Ref.  5).  In  this  case,  the 
accuracy  of  the  measurement  is  exceedingly  good,  and  is  limited  principally  by 
the  repeatability  of  the  setup,  accuracy  of  the  standards,  quantization  in  the 
sampling,  etc.  However,  the  coherent  errors  cannot  always  be  determined 
through  calibration,  and  in  some  situations,  e.g.,  multipath,  the  phasing  of 
the  error  component  may  be  time-varying.  These  statistics  are  also  useful  in 
specifying  the  component  performance  needed  to  achieve  overall  system 
performance . 


II.  STATISTICAL  ANALYSIS 


This  section  develops  the  necessary  analysis  for  the  first  and  second 
order  statistics  for  power,  voltage,  and  phase  measurements.  The  resulting 
statistics  will  be  compared  to  the  more  familiar  Gaussian  error  statistics. 

A.  STATISTICS  FOR  POWER 

The  error  statistics  for  power  measurements  are  straightforward  to 
derive.  With  reference  to  Fig.  1,  the  total  power  measured  is  given  by 

P  =  (1  +  a  eja)  (1  ♦  a  eja)»  (1) 

,  2 

=  1  +  a  +  2a  cosa 

The  mean  value  of  the  resultant  power  is  obtained  from 

2n  2 

Ep  =  (1/2ir)  J  (1  +  a  +  2a  cosct)da  (2) 

o 


The  true  power  level,  which  would  be  observed  when  a  equals  zero,  is  unity,  so 
that  the  mean  power  error  equals  a  .  Similarly,  the  variance  of  the  resultant 
power  is  obtained  from 


_2tt 


V  =  (1/2ir)  J  (P  -  Ep)2  da 


2a 


(3) 


and  the  corresponding  standard  deviation  is  given  by 

op  =  /2  a 


(4) 


T 


Typical  values  for  the  error  statistics  are  plotted  in  Fig.  2  in  the  form  of 
dB  values.  The  peak-to-peak  levels,  which  are  commonly  seen  as  the  modulation 
envelope  of  coherent  interaction,  are  presented.  The  mean  error  has  an 
appreciable  value  as  a  increases.  Because  of  the  logarithmic  form  of  the  dB 
values,  a  ±  1  o  spread  about  the  mean  value  of  the  power  is  presented  in  this 
figure. 
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Fig.  2.  Coherent  Power  Error  Statistics 
B.  STATISTICS  FOR  VOLTAGE 

The  statistics  for  voltage  errors  are  not  as  simply  derived  as  the  power 


errors.  The  total  voltage  is  the  square  root  of  Eq.  1,  and  the  error  is 
obtained  subtracting  unity,  the  true  value.  The  mean  value  of  the  total 
voltage  is 


V. 


E  =  (1/2it)  J  V  da 


=  ( 2/n )  (1  +  a)  E[4a/(1  +  a)  ] 

where  E(  )  is  the  complete  elliptic  integral  of  the  second  kind.  While  the 
elliptic  integral  is  tabulated,  a  simple  analytic  expression  is  more  useful. 
Two  polynomial  approximations  for  the  elliptic  integral  (Ref.  6)  were 
examined;  the  polynomial  series  were  derived  for  the  argument  and  the 
complement  of  the  argument.  The  series  based  on  the  argument  results  in 

Ey  -  [(1  ♦  a)  ( 1  -  (a/(  1  +  a)2  -  (3a2/4(1  -  a)V..]  (6) 

The  series  based  on  the  complementary  argument  results  in 

Ey  -  1  ♦  (a2/U)  ♦  (aH/6U)  ♦  .  .  .  (7) 

Numerical  comparisons  of  the  expressions  derived  for  the  expected  value  of  the 
total  voltage  were  made.  The  elliptic  integral  was  computed  from  an  IMSL 
subroutine  based  on  Chebyshev  polynomials  available  on  our  computer  system, 
numerical  integration  of  the  integral  using  a  trapezoidal  rule  with  1000 
points,  and  the  two  series  expansions  given  in  Eqs .  6  and  7.  The  results  of 
these  comparisons  are  presented  in  Table  I  in  terms  of  the  expected  value  of 
the  total  voltage;  the  expected  error  is  obtained  by  subtracting  1  from  these 
numbers.  The  exact  value  for  a  =  1  can  be  derived  analytically  and  is  equal 
to  4/n,  which  provides  another  check.  The  series  presented  in  Eq.  7  provides 
not  only  the  most  accurate  series  result,  but  also  the  most  convenient  form  to 
compute.  For  a  <  0.5  (-6  dB) ,  the  accuracy  of  the  mean  error  computed  from 
(1/4)  a2  is  within  1.6%. 

The  second  order  statistics  for  the  voltage  were  also  derived.  The 
expected  value  of  the  square  of  the  voltage  is  identical  with  the  expected 
value  of  the  power  given  in  Eq.  2.  The  variance  of  the  total  voltage  is  given 


Table  I.  Computational  Comparison  for  Mean  Value  of  the  Total  Voltage 


a 

IMSL 

Subroutine , 
Integration 

Eq.  (6) 

(1  +  a )  ( 1 - p] 

(1  +  aT 

Eq.  (7) 

1  +  (a2)/9 

.000 

1 .00000000 

1 .00000000 

1.00000000 

.100 

1 .00250157 

1.00909091 

1 .00250000 

.200 

1 .01002525 

1.03333333 

1 .01000000 

.300 

1 .02262952 

1 .06923077 

1 .02250000 

.400 

1 . 04041 709 

1.11 428571 

1 . 04000000 

.500 

1 .06354441 

1 . 1 6666667 

1 .06250000 

.600 

1 .09223858 

1.22500000 

1.09000000 

.700 

1.12682867 

1.28823529 

1  .12250000 

.800 

1 .16780951 

1.35555556 

1 .16000000 

.900 

1 .21 600092 

1 .42631579 

1 .20250000 

1  .000 

1 .27323928 

1.50000000 

1 .25000000 
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V  -  1  +  a2  -  {[2(1  +  a)/ir)E(W(1  +  a)2]j 


The  computation  of  the  elliptic  integral  arises  again.  The  first  two  terms  of 
Eq.  7  result  in  the  following  expression  for  the  standard  deviation 


ov  =  C(a//2)  (1  -  (3a  /1 6)  -  ...] 


The  results  of  a  numerical  comparison  between  this  expression  and  Eq.  8 
evaluated  by  the  Chebyshev  polynomial  routine  are  presented  in  Table  II.  For 
a  <  0.5  (-6  dB) ,  the  expression  in  Eq.  9  is  accurate  within  0.1%.  Values  of 
the  peak-to-peak  variations,  the  mean  error,  and  the  rms  spread  about  the  mean 
for  the  voltage  are  plotted  in  Fig.  3- 

C.  STATISTICS  FOR  PHASE 

The  evaluation  of  the  statistics  for  the  phase  is  difficult  to  derive  in 
closed  form.  Again,  referring  to  the  phasor  diagram  in  Fig.  1,  the  error  in 
the  phase  can  be  expressed  as 

e  =  tan  1  [asina/(l  +  acosa)]  (10) 

The  mean  value  of  the  phase  error  can  be  analytically  demonstrated  to  be  zero 
by  observing  the  symmetry  of  the  integrand  used  to  calculate  the  expected 
value.  The  second  order  statistics  can  be  derived  from 


V  =  ( 1 /2n)  j  (e)2  da 


This  integral  is  difficult  to  evaluate  in  closed  form.  If  the  integrand  is 
expanded  in  a  Taylor's  series,  the  following  expression  results: 

V,  -  a2/2  -  a4/8  +  ...  (1 


i 


i 


Table  II.  Computational  Comparison  for  Standard  Deviation 
of  Voltage  Errors 


a 

IMSL  Subroutine 

Eq.  (9) 

.000 

.00000000 

.00000000 

.100 

.070641123 

.07064436 

.200 

.14088642 

.14089003 

.300 

.21030662 

.21033456 

.  -400 

.27844620 

.27856777 

.500 

.34478005 

.34516753 

.600 

.40867453 

.40969501 

.700 

.46931562 

.47168912 

.800 

.52556726 

.53065997 

.900 

.57562294 

.58608073 

1  .000 

.61551690 

.63737744 

PEAK-TO-PEAK  ERROR 


Ams  ERROR  SPREAD 
(±  1  standard  deviation 
about  mean  value) 

■MEAN  ERROR  \ 


0.01  L 

-10  -20  -30  -4( 

RELATIVE  COMPONENT  ERROR  LEVEL,  dB 

Fig.  3»  Voltage  Error  Statistics 

and  the  standard  deviation  is  obtained  from  the  square  root  of  these  values. 

Table  III  shows  a  comparison  of  the  series  approximation  with  the  numerical 

evaluation  of  the  integral  for  the  standard  deviation  values.  For  a 

<  0.5  (-6  dB) ,  the  accuracy  of  the  series  form  compared  to  the  numerical 

integration  is  better  than  0.5$.  The  case,  a  =  1,  can  also  be  integrated 

2 

exactly,  with  the  result  it  /1 2  which  provides  another  check  on  the  results. 

Values  for  the  peak-to-peak  errors  and  the  standard  deviation  for  the 
phase  are  also  presented  in  Fig.  A.  The  peak-to-peak  errors  are  the 
modulation  envelope  of  the  phase  ripple  observed  experimentally. 

D.  SUMMARY  AND  COMPARISON 

For  convenience,  the  expressions  for  the  mean  and  standard  deviation  of 
the  errors  are  presented  in  Table  IV.  The  approximate  expressions  derived 


from  series  expansions  are  indicated  by 


These  statistics  are  for  the  error 


values;  the  correct  value  of  1  has  been  taken  out. 


•vv*: 


’"V*V jm-  'V-W 


m\ 


Table  III.  Computational  Comparison  for  Standard  Deviation 
of  Phase  Errors 


a 

Numerical 

Integration* 

Eq.  (12) 

0 

0 

0 

0.1 

4.08 

4.06 

0.2 

8.24 

8.14 

0.3 

12.31 

12.29 

0.4 

16.59 

16.53 

0.5 

20.98 

20.88 

0.6 

25.92 

25.38 

0.7 

30.12 

30.05 

0.8 

35.92 

34.91 

0.9 

42.54 

39.98 

1.0 

51 .96 

45.30 

* 

Values  in  degrees 

Table  IV.  Error  Statistics 

Power  Voltage  Phase 

Mean  a2  *  a2/4  0 

Standard  Deviation  /2  a  =  a//2  ■=  a//2 


Gaussian  (+5  o) 
Standard  Deviation 


0 . 4  a 


0.2  a 


0.2  a 


•  I  I _ I _ I _ I _ I _ 1 _ I 
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RELATIVE  COMPONENT  ERROR  LEVEL,  dB 


Fig.  H.  Phase  Error  Statistics 

The  comparison  between  the  statistical  values  derived  for  the  coherent 
errors  with  those  derived  from  the  Gaussian  statistics  was  done  in  the 
following  manner.  The  standard  deviations  for  Gaussian  statistics  are 
typically  derived  by  assuming  the  peak-to-peak  variations  represent  ±  5o 
values,  with  the  results  for  power,  voltage,  and  phase  given  in  Table  IV.  The 
standard  deviations  for  the  coherent  error  case  are  significantly  higher  than 
those  for  the  incoherent  case.  In  addition,  the  coherent  error  statistics  for 
both  power  and  voltage  have  non-zero  mean  errors  in  contrast  to  the  incoherent 
case . 

The  theoretical  basis  between  coherent  and  incoherent  error  statistics 
also  provides  a  contrast  between  the  two  types  of  error  components.  The 
statistics  for  incoherent  errors  are  derived  under  the  assumption  of  a  large 
number  of  statistically  similar  components  having  a  zero  mean  error.  The 
Central  Limit  Theorem  is  then  invoked  to  obtain  Gaussian  statistics 


i 

I 

4 


15 


(Ref.  7).  As  will  be  discussed,  many  practical  microwave  measurement 
applications  do  not  have  a  "large"  number  of  statistically  similar 
components.  The  question  of  "largeness"  has  been  explored  further  (Refs.  8 
and  9).  These  analyses  can  be  applied  in  the  case  of  the  statistics  for 
voltage,  and  are  based  on  the  random  walk  problem.  The  random  walk  problem  is 
typically  approached  in  a  series  of  steps,  with  the  overall  distance,  which 
corresponds  to  voltage  in  this  context,  expressed  by  Kluyver's  formula  (Ref. 
10).  Kluyver's  formula  gives  the  probability  of  the  resultant  distance  as  a 
product  of  Bessel  functions  whose  arguments  depend  on  the  individual  step 
size.  In  the  problem  depicted  in  Fig.  1,  two  steps,  1  and  a,  are  taken.  It 
can  be  shown  that  the  same  series  given  in  Eq.  7  is  obtained  from  the  random 
walk  treatment  of  the  problem. 


■ur* ■ **  «>  * 
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III.  APPLICATIONS 


Typical  applications  for  the  statistical  analyses  which  have  been 
developed  are  described  here.  This  discussion  concentrates  on  those 
measurement  applications  in  which  the  full  vector  characterization  and 
calibration  of  the  measurement  system  as  used  in  modern  network  analyzer 
systems  are  not  practical;  such  examples  are  seen  in  RF  measurements  performed 
in  open  systems  and  three  such  applications  are  described.  The  manner  in 
which  the  receiver  processes  information  also  dictates  the  appropriate 
application  of  the  statistics  for  power,  voltage,  or  phase. 

A.  RF  FACILITY  LIMITATIONS 


Antenna  and  RCS  measurements  are  performed  in  open  facilities  whose 
inherent  background  returns  limit  measurement  accuracy.  These  facilities, 
whether  outdoor  ranges  or  indoor  anechoic  chambers,  strive  to  provide  a  free 
space  environment  for  the  antenna  or  radar  target  under  test,  and  the 
measurement  errors  for  the  facility  result  from  the  deviation  from  an  ideal 
free  space  environment.  The  facility  is  illuminated  by  the  same  signal  used 
to  evaluate  the  antenna  or  radar  target,  and  consequently,  the  errors  are 
coherently  related  to  the  desired  signal. 

Reflectivity  is  the  measure  of  the  facility  background  for  antenna 
measurements .  The  physical  basis  of  this  term  is  the  ratio  of  the  reflected 
energy  from  the  facility  to  the  energy  in  the  desired  direct  signal.  The 
mechanisms  of  reflectivity  may  be  either  distributed,  as  in  the  walls  of  an 
anechoic  chamber,  or  discrete,  as  in  reflecting  objects  in  outdoor  ranges. 

The  application  of  this  term  is  illustrated  by  the  pattern  measurement 
depicted  in  Fig.  5a.  Reflectivity  is  referenced  to  the  peak  gain  level  of  the 
measured  antenna  pattern.  When  the  reflectivity  equals  the  dynamic  range  of 
the  pattern  measurement,  the  facility  reflects  as  much  power  as  that  received 
by  the  true  antenna  pattern  measured  in  an  ideal  free  space  environment. 
Generally,  antenna  pattern  accuracy  is  specified  at  an  absolute  gain  level, 
and  the  difference  between  that  level  and  the  reflectivity  level  of  the 
facility  determines  the  value  of  "a"  as  used  here.  In  the  illustration,  the 
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reflectivity  level  is  indicated  as  a  constant  level;  when  discrete  reflection 
sources  exist  in  a  facility,  the  reflectivity  varies  with  aspect  angle. 


The  RCS  facility  background  is  somewhat  more  complex.  Three  mechanisms 
contribute  to  the  background  level  of  the  facility:  (1)  the  radar  return  from 
the  physical  structure  of  the  facility,  (2)  the  radar  return  from  the  target 
support  system,  and  (3)  the  finite  isolation  between  the  instrumentation 
radar's  transmitter  and  receiver.  These  three  mechanisms  are  treated 
independently  in  error  budgets,  and  referenced  to  an  absolute  RCS  level  at  the 
range  at  which  the  target  is  measured.  RCS  measurement  accuracy  is  generally 
specified  at  an  absolute  RCS  level,  and  the  value  of  "a"  is  determined  for 
each  mechanism  relative  to  that  specified  level  as  depicted  in  Fig.  5b.  The 
aggregate  error  for  the  three  independent  mechanisms  can  be  obtained  from  an 
rss  sum  added  to  the  sum  of  the  mean  values.  This  process  can  be  easily 
justified  by  generalizing  the  power  statistics  to  N  phasors  appropriate  to  the 
definition  of  RCS.  Therefore,  this  combination  of  the  three  mechanisms  to 
form  an  aggregate  background  level  of  the  RCS  facility  is  valid. 

Both  antenna  gain  and  RCS  are  defined  on  the  basis  of  power,  and  the 
power  error  statistics  previously  derived  can  be  directly  applied.  Voltage 
and  phase  statistics  have  more  limited  application  for  RF  facility 
limitations;  one  example  lies  with  RCS  measurements  in  which  the  polarization 
properties  of  radar  targets  are  to  be  transformed  from  one  orthogonal  set  to 
another,  i.e.,  from  linear  to  circular.  In  this  case,  amplitude  and  relative 
phase  measurements  are  required,  and  a  scattering  matrix  treatment  developed 
by  Kennaugh  (Ref.  1 1  )  is  used.  The  statistics  for  voltage  and  phase  can  be 
applied  to  assess  measurement  accuracy  in  this  case. 

B.  MULTIPATH 


Statistical  analyses  have  been  widely  used  to  project  the  performance  of 
RF  systems  that  operate  in  a  multipath  environment.  A  statistical  treatment 
for  multipath  is  particularly  appropriate  because  in  many  applications  the 
phasing  between  the  multipath  components  is  time-varying.  The  time  variation 
can  result  from  changes  in  the  link  geometry,  changes  in  refraction 
characteristics,  or  in  the  case  of  ionospheric  propagation,  changes  in  the 
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plasma  profiles.  The  statistical  treatment  for  such  problems  has  used 
Rayleigh  or  Rician  distributions  (Ref.  12).  These  statistical  treatments 
provide  an  interesting  contrast  to  the  coherent  error  analyses  presented  here. 

The  distinction  between  the  coherent  error  analyses  and  Rayleigh  or 
Rician  statistics  lies  with  the  number  of  components.  As  formulated  here,  the 
statistics  consider  the  case  in  which  two  components  comprise  the  total 
field.  By  contrast,  Rayleigh  statistics  assume  an  infinite  number  of 
components  with  the  same  statistical  distribution,  while  the  Rician  statistics 
assume  one  dominant  component  together  with  a  Rayleigh  distributed  infinite 
number  of  components  (Ref.  12).  Indeed,  both  distributions  can  be  derived  by 
expanding  the  number  of  terms  in  Kluyver's  formula  used  in  the  random  walk 
problem.  More  general  distributions  are  treated  in  Ref.  13.  In  contrast  with 
these  statistical  assumptions,  many  practical  microwave  problems  have  only  a 
limited  number  of  error  components.  The  treatment  of  a  limited  number  of 
components  has  received  less  attention  (Refs.  8  and  9).  The  numerical  results 
suggest  at  least  five  or  six  components  are  required  for  the  probability  to 
approach  the  Rayleigh  distribution. 

The  statistical  analyses  for  the  coherent  errors  treat  a  direct  signal, 
represented  by  unit  amplitude,  and  one  multipath  component,  represented  by  the 
relative  amplitude  a.  The  statistical  analyses  therefore  apply  to  the  classic 
geometric  optics  picture  of  propagation  over  a  smooth,  flat  ground.  The 
Fresnel  reflection  coefficients  (Ref.  1U)  adjusted  by  any  differential  space 
loss  and  antenna  pattern  factor  can  be  used  to  estimate  the  value  of  a  for  a 
given  polarization.  Again,  the  appropriate  statistical  value  can  be  selected 
based  on  the  manner  in  which  the  receiver  processes  the  RF  signal. 

C.  ANTENNA  CROSS  POLARIZATION  ERRORS 

The  cross  polarization  inherent  in  practical  antenna  designs  leads  to 
other  applications  for  the  statistical  analyses  developed  here.  The  cross 
polarization  errors  for  a  single  antenna  are  one  example.  The  measurement  of 
the  power  and  polarization  characteristics  of  an  incident  field  by  two 
nominally  orthogonally  polarized  antennas  is  a  second  example.  Polarization 
concepts  are  described  in  more  detail  in  Ref.  15. 
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Pattern  measurements  of  antennas  typically  characterize  the  relative 
amplitudes  of  the  principal  and  cross  polarized  components.  Measurements  of 
the  relative  phase  between  principal  and  cross  polarized  components  can  be 
performed  albeit  with  additional  measurement  expense  and  calibration 
uncertainty.  The  analysis  presented  here  can  be  used  by  equating  the  relative 
level  of  the  cross  polarized  response  to  the  parameter  a  used  in  the 
statistical  analysis.  The  choice  of  power,  voltage,  and/or  phase  statistics 
again  depends  on  the  manner  in  which  the  receiver  processes  the  information. 

The  second  application  related  to  cross  polarization  for  this  statistical 
analysis  lies  with  those  systems  designed  to  measure  the  polarization  of 
incident  fields.  Ideally,  the  characterization  is  performed  with  two 
orthogonally  polarized  antennas;  in  practice,  some  level  of  cross  polarization 
exists  and  the  two  antennas  are  only  nominally  orthogonally  polarized.  In 
some  cases,  the  voltages  received  by  these  two  antennas  with  their  inherent 
cross  polarization  can  be  expressed  as 
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where  V1  and  V2  are  the  voltages  at  the  nominally  orthogonal  output  ports  of 

the  antenna,  r  is  the  complex  value  of  the  antenna  cross  polarization,  and  b 
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and  /  1  -  b  /B_are  the  incident  field  components  resolved  into  orthogonal 
components.  The  incident  field  is  specified  in  this  fashion  to  express  the 
most  general  elliptically  polarized  incident  field  on  the  two  antennas.  A 
simple  example  of  such  a  system  which  obeys  Eq.  13  is  two  orthogonal  dipoles 
connected  through  a  hybrid  network. 

If  the  two  antennas  were  ideally  orthogonally  polarized,  the  total  power 
in  the  incident  field  could  be  obtained  by  simply  summing  the  power  in  each 
antenna  port.  The  power  received  by  the  two  ports  is  given  by 

P  =  i  +  r2  ♦  Ab  /  1  -  b2  cosB  Rer  (1^ 
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It  can  be  shown  that  the  worst  case  errors  in  the  power  measurement  occur  when 
b  *=  1 //2^  and  cosg  ±1*  Under  these  conditions,  this  power  expression  reduces 
to  Eq.  1.  Physically,  these  two  conditions  correspond  to  ideally  linear  or 
circularly  polarized  incident  fields  transmitted  with  either  the  same  or 
opposite  phasing.  It  is  interesting  to  note  that  these  are  the  most  common 
choices  of  orthogonally  polarized  fields. 

The  measurement  of  an  incident  field  which  is  ideally  linearly  polarized 
is  another  typical  task.  The  polarization  characteristics  of  the  antenna,  the 
parameters  1  and  f  in  Eq.  13,  are  commonly  specified  in  either  orthogonal 
linear  or  orthogonal  circular  components.  The  incident  power  and  orientation 
of  the  electric  field,  referred  to  as  tilt  angle,  are  typically  the  quantities 
to  be  determined  by  such  systems.  The  errors  in  such  a  measurement  are 
described  in  Table  V.  The  orientation  or  tilt  angle  of  the  field  is  expressed 
by  6,  which  is  defined  0  to  180°.  The  functional  form  of  the  errors  in  Table 
V  therefore  reduce  to  the  same  form  as  those  used  in  the  statistical  analysis 
previously  described.  The  power  errors  are  identical  to  those  previously 
discussed  and  the  tilt  angle  errors  correspond  to  those  for  the  phase  error 
analysis,  and  the  analyses  presented  here  also  apply.  The  errors  for  more 
generally  polarized  incident  fields  are  described  in  Ref.  16. 

D.  VSWR  INTERACTION 

The  VSWR  interaction  between  connecting  electronic  components  is  the 
final  application  discussed  here.  Modern  network  analyzer  systems  use  vector 
calibration  to  minimize  such  errors;  however,  such  instrumentation  is  not 
always  available.  The  statistical  approach  also  finds  application  when 
component  requirements  must  be  specified  in  an  error  budget  projecting  overall 
operational  system  performance. 

A  first  order  VSWR  interaction  is  shown  in  Fig.  6.  Reference  planes 
associated  with  the  two  electronic  components  are  indicated  by  the  subscripts 
and  the  two  components  are  assumed  to  be  connected  by  a  cable  of  length  d  and 
phase  velocity  k.  The  ratio  of  the  direct  transmission  to  the  first 
interaction  as  derived  in  this  figure  is  given  by 


Table  V.  Cross  Polarization  Errors  for  Linearly  Polarized  Field 
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This  derivation  provides  both  amplitude  and  phase  information  and  the 
interaction  effects  are  deterministic.  In  many  cases;  however,  only  VSWR  data 
are  available  as  a  description  of  the  terminals  and  the  phase  information  is 
lost.  In  such  cases,  the  statistical  description  developed  here  is 
appropriate.  The  value  of  a  expressed  in  terms  of  the  scalar  VSWR  parameter 
is  given  by 

VSWR1  -  1  V SWR2  -  1 

=  VSWR1  ~  1  VSW R2  “  i  (l6) 

where  the  subscripts  on  the  VSWR  quantities  are  associated  with  the  two 
terminal  planes.  The  parameter  |fi|  can  be  equated  to  a  and  the  statistical 
analyses  can  be  applied. 

The  analysis  thus  far  has  considered  only  the  first  order  interactions; 
the  higher  order  interactions  become  significant  only  when  both  VSWR  values 
are  large.  The  effects  of  these  higher  order  interactions  are  conventionally 
written  in  closed  form  by  observing  that  each  succeeding  interaction  equals 
the  preceding  interaction  multiplied  by  R.  The  total  VSWR  interaction  is  thus 
expressed  by 

R’  =  R/(1  -  R)  (17) 
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This  expression  can  also  be  inverted  to  yield 


R  «  R'/(1  +  R' ) 
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Examining  the  phase  of  these  expressions  yields 


p  =  p'  -  p" 


where 


p"  =  tan  1  [  ( j  Ft  |  sinp)/ ( 1  -  | R  j cosp )  ] 

=  tan  1  [ ( |R' jsinp' )/(1  +  |R'|cosp')] 


and  thus  the  statistics  of  the  phase  of  the  higher  order  interactions  have  the 
same  functional  form  as  the  phase  statistics  previously  derived. 

The  effect  of  the  magnitude  of  the  higher  order  interaction  terms  can  be 
treated  statistically.  The  magnitude  of  R'  is  given  by 

| R ' |  =  1 R  I / ( 1  +  | R | 2  -  2jR|cosp)1/2  (20) 


The  highest  level  reflection  value  is  given  by 


iH'lmax  =  iR!/(1  -  lRl) 


(VSWR1  -  1)(VSWR2  -  1)/2(VSWR1  +  VSWR2) 


This  represents  the  worst  case  VSWR  interaction  and  is  plotted  in  Fig.  7  for 
representative  VSWR  values.  The  minimum  level  reflection  value  is  given  by 
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The  mean  value  of  the  VSWR  interaction  can  be  calculated  by  assuming  p  is 
uniformly  distributed  from  0  to  360°.  The  result  is 


[ 2  j  R | / ( 1  +  )R| )]K[2/[r]/(1  +  | R | ) ] 


where  K(  )  is  the  complete  elliptic  integral  of  the  first  kind.  Series 
expressions  for  the  elliptic  integral  were  examined  (Ref.  6),  and  a  form  based 
on  the  complementary  argument  was  found  both  accurate  and  computationally 
convenient;  the  expression  is 

E|R1 1  -  R[  1  ♦  ( j  R  i  2/M)  +  (9|R|V6JJ)  ♦  ...]  (24) 

The  accuracy  of  this  expression  is  better  than  I.Of  for  a  <  0.5  (-6  dB) ,  and 
also  provides  a  convenient  means  to  assess  the  effects  of  higher  order 
interactions  in  an  average  sense.  Numerical  values  for  the  effects  of  the 
higher  order  interaction  statistics  are  presented  in  Fig.  8. 
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REFLECTION  STATISTICS  WITH 


IV.  SUMMARY 


The  statistics  for  coherent  power,  voltage,  and  phase  errors  are  derived 
under  the  assumption  that  the  amplitude  of  an  error  component  has  a  constant 
level  and  that  the  phase  of  the  component  is  unknown  and  equally  likely  and 
uniformly  distributed  between  0  and  360°.  These  statistics  differ  from  the 
values  projected  from  Gaussian  statistics  which  apply  to  incoherent  errors. 

In  many  practical  microwave  measurements,  a  large  number  of  statistically 
similar  components  is  not  present,  and  the  coherent  errors  are  better 
characterized  by  this  statistical  model.  Several  applications  are  cited  in 
which  the  vector  calibration  necessary  to  determine  the  phase  angle,  a,  is 
impractical . 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -field-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  rauiatLuu  euects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 


Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 


Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  raicrowave/ratllimeter  wave  measurements, 
diagnostics  and  radiometry,  mi crowa ve/mi 1 1 imeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 


Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  ancf  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 


Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  inter  ictions ,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electro  ignetlc  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


